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Abstract 



o 

~l , We investigate the influence of nearby two-level systems on the dynamics of a qubit. The intrinsic 

P^ ■ decoherence is given by a coupling of both the qubit and the two-level systems to a heat bath. Assuming 

weak interactions between the qubit and the two-level systems we show that the dephasing of the qubit 
C^ , follows a multi-exponential decay. For a flat distribution of energy splittings of the two-level systems the 

^ ' multi-exponential behavior simplifies to a good approximation to a single exponential decay. The rate 

is given by the sum of the dephasing rate of the isolated qubit plus all relaxation rates for the two-level 
C^ ■ systems out of their present state. This latter contribution sums up to a rate approximately given by 

the mean coupling between the qubit and the two-level systems times the number of thermally active 

neighbors. 

^' 

rj . In recent years, extensive studies have been done on quantum computing. One of the major problems to 

O ' overcome in that field is the decoherence of the qubit dynamics. Besides the dephasing due to experimental 

, ^, . imperfections the interactions with the environment yields an intrinsic source of decoherence ||l|, g, |^, Q]. 

The so-called quality factor of the qubit, the number of quantum operations performed during the qubit 
coherence time, should be at least 10^ to allow for quantum error correction^. Wc will not discuss the 
^ ' various physical sources of decoherence. Instead we investigate a model where the qubit is weakly coupled to 

defects described as two-level systems (TLSs). Since the qubits are also described as TLSs, we do formally 
not distinguish between qubits and defects. All TLSs/qubits are coupled to a heat bath represented by a set 
r~^ , of harmonic oscillators providing the basic source of relaxation and decoherence. Both result from the energy 

^^ ' exchange with the heat bath like in the spin-boson model which is different from the model proposed by 

Unruh R] . Our investigation focuses on the question how neighboring TLSs influence the dephasing behavior 
of a given qubit. In the following we use qubit to address the TLS of interest and TLS for all the systems 
which form a background for the qubit. Qubits which are not used for an actual calculation contributes in 
that picture as a background. 
We restrict our investigation to longitudinally coupled TLSs leading to the Hamiltonian 

^ A 



with the Pauli matrices ai, the energy splittings A^ and the interaction constants Jij. We allow for distributed 
splittings and couplings. The many-body eigenstates are given by the product states of the isolated TLSs. 

}J] ' Nevertheless, the longitudinal coupling changes the energies of the excitations, thus leading to entanglement 

5^ \ between the quantum states of the noninteracting qubits. 

The bosons of the heat bath, Hb = 'J2k ^kb\J>k, couples with their strain field linear to the dipole operator 



ai of the two-level systems 



^SB-VaW-VA^ (6. + &!,). (2) 
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This coupling reflects the natural picture that the two-level system can change its state via an absorption 
or an emission of a boson. 



For a single qubit this is the widely investigated spin-boson model[^, 0|. The case of a transversely coupled 
pair is investigated by Terzidis et al. ||] and Dube et al. H. It is well know that a coupling to bosons leads 
to an additional transversal coupling between the two-level systems of the form Jijai" crl through virtual 
phonon exchange |12, O. We assume in the following that we can neglect that additional coupling. We 
are mainly concerned in the limit of weak coupling between the qubits and the heat bath. We focus on the 
dissipative dynamics generated by the coupling to the environement and thereby, on the dephasing times. 
The renormalization of the excitation energies are beyond the scope of the present paper. We treat the 
coupling within a resumed perturbative approach |lG, ni| as introduced by Wiirger which is a weak coupling 
approach. 

For sufficiently dilute two-level systems, the time evolution of the boson operators is not affected by the 
two-level systems and the effect of the heat bath is entirely characterized by the coupled density of states 

J,M = ^^|A«p<5(c.-c.,) (3) 



k 



with the cut-off frequency ujc of the bath and a dimensionless coupling constant a*-*' . The right hand side of 
eq. d3) assumes a linear dispersion of the bosons with an sharp upper cut-off. An exponent s = 1 describes 
systems exhibiting ohmic dissipation H and accordingly the bath is called ohmic. An exponent s > 1 is 
usually called a supcr-ohmic bath where especially s = 3 is relevant for a defect in a solid interacting with 
acoustic phonons. 

We will start with a discussion of a pair of TLSs. Its spectrum as plotted in figure 1 consists of four levels. 
Next to the levels we give the states as product states of the ground (down) and the excitated (up) state of 
the two TLSs. The full arrows give the transitions induced by fields coupling to the dipole of the first TLS 
and the dashed arrows the ones resulting from the second TLSs. 

Let's assume that TLS 1 is our qubit and we prepare it in a superposition between up and down. Accordingly 
TLS 2 is a defect. For an isolated qubit the decay of coherence is given by the probability of a boson absorption 
or emission. At temperature below the energy splitting Ai of the qubit, phonons with energy Ai are rare 
and the main effect is given by spontaneous emission. However, for a pair of coupled TLSs it is necessary to 
have information on the state of the second TLS. If it is in its down-state, boson absorption of the second 
qubit would also lead to decoherence for the first qubit since the original (coupled) state is altered. If the 
second qubit is in its up-state, boson emission of the second qubit contributes to the decoherence of the 
first one. This reflects the entanglement between the two TLSs. For two entangled qubits performing a 
quantum calculation, it is obvious that the calculation is lost if one of the TLSs dephases. Accordingly the 
dephasing rate of the pair is the sum of the rates for the two isolated qubits. (This argument holds only 
true as long as the interaction does not change the quantum mechanical evolution of the qubits meaning 
J <C A;.) Since the above argument was independent of the strength of the interaction between the TLSs, 
the question arises, how we can avoid entanglement of a qubit with surrounding TLSs which do not take 
part in the actual computation. 

To answer this question we derive the memory kernel, resulting from the coupling to the heat bath, for the 
dynamics of the dipole operator ai of the first qubit. Following the approach in ref. |1G] we calculate the 
memory kernel within the lowest contributing order of the coupling between system and bath and apply a 
Markov approximation. Neglecting the real part of the memory kernel, which represents the renormalization 
of the energy splittings, we end up with the rates for the various dynamic processes. We restrict ourselves 
to the subspace of the operators, |ti)(ii| and |tT)(iTI; which correspond to a transition from down to up 
for the first qubit or, in other words, to a preparation in a superposition between up- and down-state. The 
dynamics of this subspace is independent from the rest space and the time evolution of the subspace of the 
complex conjugate operators is identical. 

The dynamics of the operators | ti)(ii | and | tT)(iT I are not independent from each other since the off- 
diagonal components of the memory kernel are non zero. 
The time-evolution operator in Laplace space for this subspace (in the basis of the above given operators) is 



given by 

U{z) = (4) 

'.-u+ + i{Ti + 2T2]) -J2r2x 

~i2V2] z-ijj^+i{Ti + 2T2i) 

with the excitation energies w^ = Ai =p 2 J « Ai, the usual onc-phonon rate of the first qubit Fi = 

Ji(Ai) coth(/?Ai/2), the decay rate of the second qubit due to phonon emission T21 — J2(A2)(H-n(A2)) with 

the Bose factor n{/S.2) and the decay rate of the second qubit due to phonon absorption r2| = J2(A2)n(A2). 

Thereby Ji/2(Ai/2) are the coupled density of states for the two TLSs (compare eq.(H)). The phase coherence 

rates are given by the imaginary parts of the two poles of the time evolution operator. 

This dynamics has two simple extremes. When the difference of the two excitation energies |cj+ — a;_| = 4 J is 

bigger than the off-diagonal components, we can neglect the off-diagonal entries. In this case, the dynamics 

of qubit 1 has two different rates which depend on the state of qubit 2. Therefore, the rates are determined 

by both processes: decay of the qubit 1 and decay of the qubit 2. 

In the second case, J < max{r2x,r2|} = r2x, we have to taken into account the off-diagonal entries. The 

off-diagonal elements cause that the eigenstates of the time evolution do not correspond to the eigenstates of 

the Hamiltonian. Instead the dynamics of the two qubits completely decouples and the coupling is negligible 

(for the dynamics of the qubit 1) as long as it is smaller than the linewidth of the states of the isolated qubit 

2 which is given by (r2| -I- r24)/2. 

Thus, only when the coupling between a qubit and a second TLS is smaller than the linewidth of the second 

TLS the dynamics of the qubit is independent from the second TLS. We should point out that the actual 

quantum mechanical evolution of the states are only weakly affected since J/ A ^ 1. But the dephasing rate 

might be changed considerably since we have no knowledge about ri<r2|/|. 

In the case of many qubits the above conclusion remains qualitatively valid. The decoherence rate of a qubit 

of interest depends on the actual states of all other TLSs. If we have a distribution of interaction energies, 

we have to devide all TLSs into two groups depending on the condition Jai < r^ or Jai > Ti^ where a is 

the index of the qubit of interest. 

If Jai < Til ^^^ ^~^^ '^-'^^ does not influence the qubit a. Nevertheless a flip of the qubit i will change the 

resonance energy of the qubit a. Since that change is smaller than the linewidth, we would not be able to 

measure it. In the case of many neighbors fulfiling that condition, the changes might add up resulting in a 

change in the resonance energy which is bigger than the linewidth. In an echo experiment at an ensemble of 

qubits, this leads to an 'effective' dephasing where the relative phase between the qubits decays. This effect 

is called 'spectral diffusion' Q. 

All neighbors in the second group with Jai > Tii form an entangled cluster with the qubit a. As in the 

case of a pair, the decoherence rate of a superposition of eigenstates of the qubit a, with all the TLSs in an 

eigenstate, gets: 

-p _ -p „ Y^ / rjT TLS i in the down-state ._. 

* ~ " + ^ I r,x for TLS i in the up-state ^ ' 

{ilJai>rii} 

Thereby, we assumed a clear separation between Jia and F^j^. If both arc of the same order of magnitude, 
the time evolution is a complicated mixture of both extremes. This assumption should at least be valid 
for the majority of the TLSs. If the TLSs are qubits like the qubit a, the decoherence rate scales with the 
number of qubits. This seems to be natural for a situation where N entangled qubits perform one calculation. 
If only one of these qubits dephases, the calculation is lost. The main point lies in the fact that they are 
entangled as long as the interaction between two qubits exceeds their linewidth. For example, if one wants to 
send quantum information from one to another quantum computer, the interaction must be bigger than the 
linewidth since otherwise, the quantum computers would dephase even before they are entangled (starting 
from an isolated initial condition) . After the transmission both computers must be separated again fulfilling 
the above condition. 

Within the picture of one qubit surrounded by TLSs, we have to specify the distributions of couplings P{J) 
and of the resonance energies Q(A) in order to get quantitative results. 



Assuming homogeneously distributed TLSs and a dipole-dipole like interaction between TLSs we obtain for 
the distribution of couplings 

PiJ) = ^ (6) 

between a minimal and a maximal value of the coupling with equal probability of both signs. J ~ NJmin 
normalizes the distribution to the number of TLSs N. Since J « (| J|) — J l'D-{Jmax / Jmin) , we refer to J in 
the following as the mean coupling. 
We want to consider a case in which all neighbors are equivalent 

0(A) = S{A - A) . (7) 

Without further knowledge over the neighbors we expect them in thermal equilibrium with the environment. 
In that case the average decay of the qubit a follows the form ( J^ ^ A) 

(e-") c. e-^"' ■ [p-e-^'' + p+ e.-^^'^ '''■''] (8) 

with the occupation probabilities of the ground/excited state p* — exp(±/3A/2)/cosh(/3A/2). 
This simple result has few remarkable features. First of all, we have no longer a simple exponential decay 
but a sum of two. The original decay of the qubit a is still present in the prefactor. The decay of the TLSs 
lead to the term in the brackets. This term does not depend on the actual decay rates of the TLSs but only 
on their mean coupling J to the qubit a and the probabilities of these neighbors to be in their ground or 
excited state. If the temperature is higher than the energy splitting of the TLSs (fc^T 3> A) both rates are 
equal to a good approximation. In the other limit {ksT <C A) only one decay is relevant and its contribution 
from the neighboring TLSs is vanishing. Thus only thermally active TLSs are affecting the dephasing of the 
qubit. Actually two decays are only seen in the case where the temperature matches the energy splitting. 
In the second step, you might consider neighbors with broadly distributed energy splittings. If we assume, 
for example, just a flat distribution (5(A) = (5o0(A„ii„ ^ A)0(A — Amaj;) between a minimal and a maximal 
value, we get 

i 

where i runs over all neighbors and Ji ~ QaJmin with the normalization factor Qo- In order to get simple 
analytical expressions we approximate the exponential by exp(— /3Ai) ~ Q{T — A^). This approximation 
fails for all systems with A^ ~ T. If the majority of TLSs, as in the case of a flat distribution of energy 
splittings, does not fulfill that condition, we expect that the approximation works pretty well and we get 

(e-") ~ e-ro*-2.7(T)t (^0) 

with J{T) = QoiksT — /S.min)N Jmin- This shows an additional dephasing rate which is approximately 

linear in temperature as long as the temperature is between the extremal values of the distribution (5(A). 

The dephasing due to nearby local defects is given by the mean coupling to the thermal active neighbors. In 

fact, the above used approximation divides the TLSs sharply in thermal active and inactive ones. Thereby 

we loose the multi-exponential behavior of the correct solution. 

We investigated the dephasing in a model of weakly coupled TLSs. The coupling was chosen only to influence 

the energies of the excitations and thus to ensure entanglement. It does not change the states in a sense that 

the eigenstates of a coupled cluster are the product states of the isolated TLSs. In order to cause dephasing 

we coupled all TLSs to a heat bath in a way that dephasing as well as relaxation is caused through energy 

exchange with the bath. 

We find that the dephasing of a qubit weakly coupled to N TLSs is no longer a simple exponential decay but a 

complicated multi-exponential decay. We investigated two simple cases. In both, we assumed a homogeneous 

distribution of TLSs and a dipole type coupling between the qubit and the TLSs. 

In the case of identical TLSs, we get two exponential decays instead of one. If the temperature is higher than 

the energy splitting of the TLSs (fc^T ':^ A), both rates are equal. In the other limit (fc^T <^ A) only one 



decay is relevant and its contribution from the neighboring TLSs is vanishing. Thus only thermally active 
TLSs are affecting the dephasing of the qubit. Actually the two decays are only seen in the case where the 
temperature matches the energy splitting. 

On the other hand, a flat distribution of energy splittings for the TLSs leads to a multi exponential decay. 
Nevertheless, we end up with one major dephasing rate where all thermally active TLSs are contributing 
to. The dephasing through local defects is thereby given by the mean interaction with all thermally active 
defects. Thus, the additional dephasing is independent from the actual decay mechanisme of the TLSs. 
The above results are valid as long as the linewidth of the TLSs is smaller then their coupling to the qubit. 
This fact opens another way to avoid dephasing due to nearby TLSs by increasing the linewidth of the 
thermal active TLSs. 
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Fig. 1: The spectrum of two coupled TLSs. The states are labeled by their composition of the up and down 
states of the two isolated TLSs. The possible phonon transitions are shown by the arrows. The full (dashed) 
arrows are transitions induced by interactions with the TLS 1 (TLS 2). 



